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ABSTRACT 

We present HST WFPC2 V and I images and GHRS UV spectrophotometry of the spectral 
regions around Lya and 01 A1302 of the extremely metal-deficient ( Z ~ Z0/4I ) blue compact 
dwarf (BCD) galaxy SBS 0335-052. All the star formation in the BCD occurs in six super- 
star clusters (SSC) with ages < 3-4 Myr. As there is no evident sign of tidal interaction, star 
formation in the BCD is probably triggered by stochastic cloud collisions in the HI envelope. 
Dust is clearly present and mixed spatially with the SSCs. There is a supershell of radius ~ 
380 pc, delineating a large supernova cavity. The instantaneous star formation rate is ~ 0.4 
M,Q yr"^ . 

Strong narrow Lya emission is not observed. Rather there is low intensity broad (FWZl 
= 20 A) Lya emission superposed on even broader Lya absorption by the HI envelope. This 
broad low-intensity emission is caused by resonant scattering of Lya photons. The absence of 
strong Lya emission may be due partly to dust absorption, but mainly to multiple scattering 
which removes Lya photons from the small HST aperture. As the HI cloud is seen nearly 
edge-on, geometrical effects may play also a role as photons escape more easily in a direction 
perpendicular to the plane than along it. 

The BCD appears to be a young galaxy, undergoing its very first burst of star formation. 
This conclusion is based on the following evidence: 1) the underlying extended low-surface- 
brightness component is very irregular and filamentary, suggesting that a significant part of the 
emission comes from ionized gas; 2) it has very blue colors (-0.34 < (V-l)o < 0.16), consistent 
with gaseous emission colors; 3) the 01 A1302 line is not detected in absorption in the GHRS 
spectrum, setting an upper hmit for N(0)/N(H) in the HI envelope of the BCD of more than 
3000 times smaller than the value in Orion. 

Subject headings: galaxies - young: interstellar matter - nebulae: HII regions 
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1. INTRODUCTION 

The formation of galaxies is one of the most fundamental problems in astrophysics, and 
much effort has gone into the search for primeval galaxies (PG). A possible definition of a 
primeval galaxy is a young system undergoing its first major burst of star formation. It is 
now widely believed that the vast majority of galaxies underwent such a phase at redshifts r^2 
or greater. In most galaxy formation scenarios, young galaxies are predicted to show strong 
Lyo; emission, associated with the cooling of the primordial gas and the subsequent formation 
of a large number of massive ionizing stars (Partridge & Peebles 1967; Meier 1976; Baron & 
White 1987; Chariot & Fall 1993). Yet, despite intensive searches, the predicted widespread 
population of Lya primeval galaxies has remained elusive (Pritchett 1994). The reasons for 
this non-detection are not yet clear. 

Several objects have been put forward as possible PG candidates, ranging from high-redshift 
radio galaxies to Lya emitters found around quasars and damped Lya systems, mainly on 
the basis of very high luminosity and star formation activity (Fontana et al. 1996; Pettini, 
Lipman & Hunstcad 1995; Yee et al. 1996). However, most of these candidate PGs already 
contain a substantial amount of heavy elements, implying previous star formation and metal- 
enrichment and hence not satisfying the above definition of PG. For example, the spectrum 
of the PG candidate discussed by Yee et al. (1996) shows a strong P-Cygni profile in CIV 
A1550, indicating the presence of heavy elements since low-metallicity stars do not have strong 
winds and show no or weak P-Cygni profiles. Moreover, even if PG candidates at high-redshift 
are discovered, it is difficult to study them in detail because of their extreme faintness and 
very compact angular size. Here we take a different approach to the PG problem. Instead of 
searching for very high-redshift galaxies in the process of forming, we look for nearby galaxies 
undergoing their first burst of star formation, and hence satisfying the above definition of a 
PG. The best candidates for such a search are blue compact dwarf galaxies (BCD). 

BCDs are extragalactic objects with M5>-18 where intense star formation is presently 
occuring, as evidenced by their blue UBV colors, and their optical spectra which show strong 
narrow emission lines superposed on a stellar continuum which is rising toward the blue, similar 
to spectra of HII regions. Star formation in BCDs cannot be continuous but must proceed by 
bursts because of several observational constraints: 1) Gas is transformed into stars at the rate 
of approximately 1 M0yr~^, so that the current burst cannot last more than about 10® yr before 
depleting the gas supply; 2) Optical-infrared colors of BCDs give ages of about 10^ yr; and 3) 
Population synthesis of UV spectra of BCDs give invariably jumps in the stellar luminosity 
function, indicative of starbursts (see Thuan 1991 for a review). 

Ever since their discovery, the question has arisen whether BCDs are truly young systems 
where star formation is occuring for the first time, or old galaxies with an old underlying 
stellar population on which the current starburst is superposed (Searle, Sargent & Bagnuolo 
1973). Thuan (1983) carried out a near-infrared JHK survey of BCDs and concluded that all 
the objects in his sample possessed an old underlying stellar population of K and M giants. 
That result was not unambiguous as the JHK observations were centered on the star-forming 
regions and the near-infrared emission could be contaminated by light from young supergiant 
stars. The advent of CCD detectors allowed to look for the low-surface-brightness underlying 
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component directly. Loose & Thuan (1985) undertook a CCD imaging survey of a large BCD 
sample and found that nearly all galaxies (> 95%) in their sample show an underlying extended 
low-surface-brightness component, on which are superposed the high-surface-brightness star- 
forming regions. Subsequent CCD surveys of BCDs have confirmed this initial result (Kunth, 
Maurogordato & Vigroux 1988, Papaderos et al. 1996). Thus, most BCDs are not young 
galaxies. In Kunth et al. (1988)' sample, the BCD I Zw 18 does not show an underlying 
component. This BCD with a metallicity of only Zq/50 is the most metal-deficient galaxy 
known (Searle & Sargent 1972). It has been discussed as a possible young galaxy, currently 
undergoing its first burst of star formation (Kunth & Sargent 1986). Hunter & Thronson 
(1995) obtained HST images of I Zw 18 and concluded that the colors of the underlying diffuse 
component are consistent with those from a sea of unresolved B or early A stars, with no 
evidence for stars older than ~ 10'' yr. 

For more than 20 years, I Zw 18 stood in a class by itself. The BCD metallicity distribution 
ranges from ~Z0/3 to ~Z0/5O, peaking at ~Z0/1O, and dropping off sharply for Z<Z0/1O 
(Kunth & Sargent 1986). Intensive searches have been carried out to look for low- metallicity 
BCDs but they have met until recently with limited success. For example, the majority of the 
BCDs in the Salzer (1989) and Terlevich et al. (1991) surveys have metallicities larger than 
Z0/IO. Several years ago, a new BCD sample has been assembled by Izotov et al. (1992, 1993) 
from objective prism survey plates obtained with the Im Schmidt telescope at the Byurakan 
Observatory of the Armenian Academy of Sciences during the Second Byurakan Survey (SBS, 
Markarian, Lipovetsky & Stepanian 1983). The most interesting feature of the SBS is its 
metallicity distribution (Izotov et al. 1992): it contains significantly more low- metallicity BCDs 
than previous surveys. It has uncovered about a dozen BCDs with Z<Z0/15, more than 
doubling the number of such known low-metallicity BCDs and filling in the metallicity gap 
between I Zw 18 and previously known BCDs. This low-metallicity sample has been used to 
determine the primordial helium abundance (Izotov, Thuan & Lipovetsky 1994, 1996), heavy 
element abundances in metal-deficient environments (Thuan, Izotov & Lipovetsky 1995) and 
study the BCD large-scale spatial distribution (Pustilnik et al. 1995). 

Here we focus our attention on the most metal-deficient BCD in the SBS sample, SBS 
0335-052. Using the Hubble Space Telescope (HST), we have obtained optical imaging and UV 
spectroscopy for the galaxy. We discuss the HST observations and their data reduction in §2. 
We argue for the youth of SBS 0335-052 in §3. In §4 we show how other observational data also 
support the youth hypothesis, and discuss the Lya problem in the BCD and its implications 
on Lya searches of high-redshift primeval galaxies. We also compare the properties of SBS 
0335-052 with those of other nearby young galaxy candidates. We summarize our conclusions 
in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

The BCD SBS 0335-052 was discovered by Izotov et al. (1990) to have an extremely low 
metallicity. With Z~Z0/41 (Melnick, Heydari-Malayeri & Leisy 1992), it is the second most- 
metal deficient BCD known, only behind I Zw 18. Its integral characteristics are summarized 
in Table 1. At the adopted distance^of 54.3 Mpc, 1" corresponds to a linear size of 263 pc. To 

^We shall adopt throughout a Hubble constant of 75 km s^^Mpc^^. 
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understand better the nature of this extraordinary object, we have obtained images and high 
resolution UV spectra of SBS 0335-052 with the refurbished HST , in the course of a larger 
study of extremely metal-deficient BCDs . The superior angular resolution of HST is ideal for 
studying this extremely compact (d25~12") object. 

2.1. Imaging 

We obtained images of SBS 0335-052 on 1995 January 27 during cycle 4, after the refurbish- 
ment mission, with the HST Wide Field and Planetary Camera 2 (WFPC2) in filters F569W 
and F791W, which we will refer to as V and I throughout the paper. Two exposures of equal 
duration were obtained in each filter to permit identification and removal of cosmic rays. The 
total exposure time was 1800 s in V and 4400 s in I. The scale of the WFPC2 is 0.102" per 
pixel. 

The data reduction followed the procedures described in Thuan, Izotov & Lipovetsky 
(1996a). Preliminary processing of the raw images including corrections for flat-fielding was 
done at the Space Telescope Science Institute (STScI) through the standard pipeline. Subse- 
quent reductions were carried out at the University of Virginia using IRAF and STSDAS 5 
Point sources in the processed images have a FWHM of ~ 2 pixels or 0.2". Cosmic rays were 
removed and the images in each filter were combined. We found that all exposures in a given fil- 
ter coregistered to better than ~ 0.2 pixels. The transformation of instrumental magnitudes to 
the Johnson-Cousins UBVRI photometric system as defined by Landolt (1992) was performed 
according to the prescriptions of Holtzman et al. (1995). The resulting measured brightness of 
the sky background is 22.6 mag arcsec"^ in V and 21.8 mag arcsec"^ in I. 

2.2. UV spectroscopy 

Spectroscopic observations were obtained on 1995 January 3-4 using the Goddard High 
Resolution Spectrograph (GHRS) on board HST, in two separate exposures, one with the 
G160M grating covering the 1210-1250A region around the Lya line of hydrogen and the other 
with the same grating covering the 1300-1350A region around the resonance line of oxygen at 
1302 A. The exposure time was 7180 seconds for each region. The 2"x2" rectangular Large 
Science Aperture was used, centered on the brightest knot (corresponding to cluster 1 in Table 
2 to be discussed later). The flux and wavelength calibrations were done at STScI through the 
standard pipeline. The spectral resolution is 0.12 A. 

3. A PROBABLE YOUNG GALAXY 

3.1. Young super-star clusters 

■^IRAF: the Image Reduction and Analysis Facility is distributed by the National Optical Astronomical 
Observatories which are operated by the Association of Universities for Research in Astronomy, Inc. (AURA), 
under cooperative agreement with the National Science Foundation (NSF). 

''STSDAS: the Space Telescope Science Data Analysis System. 
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We display in Figure la the V image with the contrast level adjusted so as to show the 
details of the high-surface-brightness star-forming regions. The latter are seen to contain 6 
bright stellar clusters all within a region of ~ 2 arcsec diameter, or ~ 520 pc in hnear size. It is 
clear from the contour map (Figure lb) that the clusters (labeled 1 to 6) are unresolved, with 
radii <25 pc. Aperture photometry of the point-like objects was performed using the package 
APPHOT in IRAF. The magnitudes were obtained using a circular aperture of radius 2 pixels, 
which is about the FWHM of the stellar profile. To determine the background, we use an 
annulus centered on the cluster covering 3-5 pixels. The precision of the photometry is good 
despite the high background because of the relative brightness of the clusters (V<22 mag). 

Table 2 shows the Vq and Iq magnitudes, and the (V - I)o colors for all six clusters. All 
magnitudes and colors have been corrected for Galactic extinction which is small in the direction 
of SBS 0335-052 (Ay =0.10 mag, A7=0.06 mag, Burstein & Heiles 1982). The magnitudes of 
the stellar clusters vary between Vq = 19 mag and Vq = 22 mag, which correspond to absolute 
magnitudes My^ between -11.7 mag and -14.7 mag. The (V - I)o colors are very blue for all 
clusters, ranging between -0.61 and 0.31. In fact, the (V - I)o color for the bluest and brightest 
cluster (No.l in Table 2) is bluer than that of the hottest main sequence O star ( -0.33 for a 
OV star with solar metallicity, Bessell 1990). Models of young (age <3-4 Myr) stellar clusters 
(Leitherer & Heckman 1995) predict (V - I)o in the range 0.0 - 0.1 for a Salpeter or a Miller- 
Scalo-type initial mass functions (IMF), with an upper mass limit of 100 M0 and a metaUicity 
of one tenth that of the Sun, in the two limiting cases of an instantaneous burst and a constant 
star- formation rate over a time interval of 3xl0®yr. Thus cluster 1 is ~ 0.7 mag too blue in 
(V-I)o as compared to the model predictions. Several effects can be thought of which make (V 

- 1)0 bluer. The most important effect is the contamination by gaseous emission, especially in 
the V band where several strong emission lines are present. To estimate this contamination , we 
have used the spectroscopic observations of Izotov ct al. (1996) who have measured for cluster 
1 equivalent widths of 178 A, 199 A, 616 A and 1138 A respectively for the H/?, [OIII]A4959, 
5007 and Ha lines. By convolving the line fluxes with the transmission curve for the V filter 
(HST WFPC2 Handbook 1995), we estimate the contribution of gaseous emission to the V 
band to be ~ 0.5 mag. The remaining small ~ 0.2 mag difference can be plausibly attributed 
to contamination by gaseous continuous free-free and free-bound emission and by a metallicity 
difference between models and observations. Stellar clusters in SBS 0335-052 have Z<Z0/41 
(stars have a lower metallicity than the gas) while Leitherer & Heckman (1995)' models were 
computed for Z = Zq/IO. A lower metallicity would result in a lower line blanketing, and 
hence to bluer colors. Thus we believe that the (V-I)o color of cluster 1 is consistent with the 
model calculations of Leitherer & Heckman (1995): the cluster was created in a starburst with 
a normal IMF (x = 1.35 where x is defined by dN/d(logM) oc M"^, M„ = IOOMq, M^^IMq) 
less than ~ 3-4 Myr ago. 

The clusters listed in Table 2 and labeled in Figure lb all lie roughly in a SE-NW direction, 
with cluster 1 at the SE edge. They all have redder (V - I)o colors than cluster 1. Figure 2 
shows that there is a systematic increase in reddening (although the decrease in brightness is 
not monotonic) of the clusters away from the brightest cluster (No. 1). The color gradient A(V 

- 1)0/ Ar is nearly constant and equal to ~ 0.47 mag arcsec^^ or ~ 0.18 mag/100 pc. Thus 
cluster 6 which is furthest away from cluster 1, at ~2 arcsec or 526 pc, is ~0.9 mag redder. We 



7 



argue that the reddening is not caused by evolutionary effects, but rather hj internal extinction 
by dust. The presence of dust patches is clearly seen in Figure 3 where we show the (V - 1) color 
map of SBS 0335-052. Black pixels denote blue colors while white pixels indicate red colors. 
Comparison with Figure la shows that cluster 1 lies outside the dust patches, but that the 
other clusters are in the same regions as the dust and are subjected to internal extinction. The 
dust patches are seen to systematically increase in importance away from cluster 1, confirming 
visually the effect found in Figure 2. There is quantitative evidence that internal extinction 
does increase with increasing distance from cluster 1. Melnick et al. (1992) have obtained with 
the ESQ New Technology Telescope (NTT) spectra of two knots in SBS 0335-052, the first knot 
containing clusters 1 and 2 and the second knot clusters 4 and 5. Although their seeing was 
excellent for ground-based observations (0.58" - 0.85" FWHM), thanks to the adaptive optics 
of the NTT, Melnick et al. (1992) did not have the necessary angular resolution to resolve 
those two knots into individual clusters as done here with HST. Those authors found that the 
extinction coefficient C(H/?) is only 0.07 dex for the first knot while it is larger for the second 
one, C(H/3) = 0.56 dex. This corresponds to a reddening difference of E(B - V) ~ 0.3 mag, 
which is just about the amount needed to account for the color difference A(V - I)o = 0.45 mag 
between clusters 1 and 4 (Table 2). Thus there is good evidence that most of the color gradient 
depicted in Figure 2 is not due to evolutionary effects but to increasing extinction away from 
cluster 1, from Ay — for cluster 1 to Ay ~ 1.8 mag for cluster 6. 

We now discuss the nature of the blue clusters in SBS 0335-052. After correction for 
internal extinction, gaseous emission and metallicity effects so that each cluster has the same 
intrinsic (V - I)o ~ 0.1 as cluster 1, then the cluster absolute magnitudes are in the range 
-14.1 < My < -11.9. This is precisely the range of absolute magnitudes found by O'Connell, 
Gallagher & Hunter (1994) for objects called "super-star clusters" (SSC). SSCs are defined by 
their combination of very small sizes and high luminosities. They have diameters less than 50 
pc, consistent with the radius upper limit of 25 pc set for our stellar clusters. O'Connell et 
al. (1994) could resolve some of their objects with the HST Planetary Camera and found that 
their half-light radii R0.5 < 3.5 pc, comparable to radii of normal star clusters and compact HH 
regions. In spite of their compactness, the luminosities of the SSCs are comparable to those 
of the physically larger super associations (My < -11). This implies that SSCs have surface 
brightnesses > 100 times those of clusters and associations in normal giant HH regions. They 
represent particularly intense periods of star formation where gas is transformed into stars at 
the rate of ~ IM0 yr^^ (Thuan 1991). From the H/3 luminosity in knot 1 (clusters 1 and 2) 
of SBS 0335-052 measured by Melnick, Heydari-Malayeri & Leisy (1992), we found that ~ 10=^ 
- 10^ massive O stars are present in that knot. Thus it appears that some BCDs make stars 
by forming SSCs. This SSC formation appears to be related to the starburst mode of star 
formation, i.e., to a sharp increase in the star formation rate of a galaxy during a relatively 
short (< 10^ yr) period. 

The advent of HST and its superior angular resolution has permitted the discovery of SSCs 
with an ever greater frequency, all of them in galaxies undergoing starbursts (see the review 
by Ho 1996). Examples include the amorphous galaxy NGC 1140 (Hunter et al. 1994), the 
BCD He2-10 (Conti & Vacca 1994), the pecuhar galaxy NGC 1275 (Hohzman et al. 1992), 
the merging systems NGC 7252 (Whitmore et al. 1993) and NGC 4038/4039 (Whitmore & 
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Schweizer 1995). Meurer et al. (1995) have studied nine starburst galaxies ranging from BCDs 
to ultraluminous merging far-infrared galaxies and also conclude that cluster formation is an 
important mode of star formation in starburst galaxies. Even in galaxies where star formation 
activity is not dominant, Barth et al. (1995) found SSCs within the circumnuclear starburst 
rings around the Seyfcrt galaxies NGC 1097 and NGC 6951. 

It has been suggested (e.g. Meurer et al. 1995), mainly on the basis of sizes and luminosities, 
that these blue SSCs will fade out to become globular clusters, once the starburst is over. 
We discuss this issue using our own HST data on stellar clusters in several BCDs in various 
evolutionary stages. Besides SBS 0335-052 which we shall argue, appears to be a young BCD 
undergoing star formation for the first time, we have also used HST to study the BCD Markarian 
996 which is at a much later stage of evolution, having undergone several bursts of star formation 
before the present one (Thuan, Izotov & Lipovetsky 1996a). Projected against and around 
Mrk 996, are many pointlike sources which can be interpreted not as young but old (> 10^° 
yr) globular clusters, as evidenced by their red color ((V-I)o ~ 0.9), absolute magnitude range 
(-10.4 < My < -6.8) and unresolved diameters (< 20 pc). In the -9.0 < My <-7.5 range, the 
cluster luminosity function in Mrk 996 is remarkably similar to that of the Milky Way, with 
slightly more clusters in Mrk 996 relative to the Galaxy in the magnitude range -10.5 < My 
<-9.0. The presence of red but not blue clusters around Mrk 996 implies that SSC formation 
does not always accompany a starburst in a BCD. If SSCs are indeed progenitors of globular 
clusters, SSC formation has occurred not with the present burst but with past ones in Mrk 
996. Furthermore, the spatial distributions of clusters in the 2 BCDs are very different. In the 
case of SBS 0335-052, the SSCs are within the star-forming regions, while in the case of Mrk 
996 the red globular clusters surround the galaxy, just as in the Milky Way. 

Since SSC formation is associated with merging events, we have investigated the environ- 
ment of SBS 0335-052. A VLA HI map of the BCD (Thuan et al. 1996b) shows it to be 
embedded in a very large flattened HI cloud, 4.3' EW by 1.5' NS in angular size or 64 kpc by 
24 kpc in linear size. There arc two peaks of HI column density, one of which is associated 
with SBS 0335-052 and the other with a fainter dwarf irregular galaxy (Pustilnik et al. 1996). 
There is a nearby companion, the face-on Scd galaxy NGC 1376 ( m^ = 12.8 mag, 2.0'xl.7' ) 
at ~9.5' W [ a(1950) = 03''34"^37.2^ 5(1950)= -05°12' 29 "] and with a heliocentric velocity of 
4159 km s~^, only 116 km s^^ larger than the velocity of SBS 0335-052. The VLA map shows 
no evident sign of tidal interaction between the HI cloud and NGC 1376. There are three other 
galaxies with angular size ~ 1' within an angular distance of 6.5' from the BCD, but they are 
all background galaxies : MPWIG at 4.3' NE with v= 10779 km s'^, A447 at 5.5' SE with v= 
33660 km s^^ and IRAS 03348-058 at 6.2' NW with v= 5556 km s-\ If NGC 1376 moves at ~ 
100 km s^^ relative to the BCD, it will take about 10^ yr to go the 150 kpc distance separating 
the two galaxies, much larger than the < 3-4 Myr age of the SSCs. Thus the formation of SSCs 
in SBS 0335-052 is probably not triggered by tidal interactions, but by stochastic collisions of 
small HI gas clouds within the large HI envelope. 

3.2. A very blue underlying component 

Figure la shows that there is a more extended lower surface-brightness diffuse component 
underlying the SSCs. This is seen better in figure 4 where the contrast has especially been 
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adjusted to display low-surface-brightness features. The underlying component extends in the 
same SEl-NW direction as the SSCs and is roughly elliptical in shape, with the SSCs located 
at the SE end. The nature of this underlying component is crucial for determining the age 
of the galaxy, whether it is undergoing its very first burst of star formation, in which case it 
satisfies the definition of a primeval galaxy, or whether the present burst is just one after many, 
in which case the underlying component contains an older stellar population. 

We have derived surface brightness and color profiles for the underlying component. Surface 
photometry was done by fitting ellipses to the isophotal contours using the task ELLIPSE in 
STSDAS^. We show respectively in Figures 5 a,b and c the V and I surface brightness and 
(V-I) color distributions as a function of log r. Here r is the equivalent radius of the best 
fitting elliptical isophote defined as r = (ab)^/^, where a and b are respectively the semi-major 
and semi-minor axes of the ellipse. In Figures 5a and 5b, it can be seen that for r > 1", the 
surface brightness distribution is well fitted by a power law (a straight fine in this coordinate 
system) of slope n ~ 2.3 ( where n is defined as I oc r~") for both V and I bands. The large 
statistical error bars in the radius range -0.45 < log r < -0.05 or OA" < r < 0.9", is due to light 
contamination from the shell structure seen clearly in Figure 4. This structure is probably the 
result of supernova explosions in the SSCs which have carved out a nearly spherical cavity in 
the ambient interstellar medium (ISM) of the BCD. As for the non-monotonically increasing 
error bars at radii log r ~ 0.3 or r ~ 2" and log r~0.7orr~5", they are probably due to light 
contamination from the SSCs at the SE side of the galaxy. Most remarkable is the (V-I) color 
profile (Figure 5c) which is nearly flat and shows extraordinary blue colors for the underlying 
low-surface-brightness component. The (V-I) colors vary between -0.3 and 0.2, corresponding 
to foreground extinction corrected colors (V-I)o between -0.34 and 0.16. 

If the underlying emission is of stellar origin, this range of colors would correspond to spectral 
types between OV and A5V stars, using the color-spectral type calibration of Bessell (1990). 
This calibration holds for solar-metallicity stars. For stars with < 1/41 of the Sun's metallicity, 
as is the case for SBS 0335-052, the corresponding spectral types at a given (V-I) color are 
slightly later. In any case, if stars are responsible for the underlying extended emission, they 
cannot be older than ~ 10^ yr. However we do not believe the underlying emission to be of stellar 
origin. The irregular, blotchy and filamentary structure of the underlying component argues for 
a gaseous origin. This hypothesis is further supported by the colors of the shell structure which 
are themselves very blue, as seen clearly in the color map in Figure 3. Aperture measurements 
at various locations along the shell give (V-I)o between -0.75=b0.20 and -0.18±0.20. The bluest 
colors are bluer than that of the hottest star ((V-I)o=-0.33 for an OV star according to Bessell 
(1990)), but can be understood if due to gaseous emission. Another line of evidence is provided 
by the long-slit observations of Izotov et al. (1996) of the underlying component. Those authors 
detected strong [OIII]A4363 emission indicating high excitation gas as far as ~ 2 kpc from the 
SSCs. Using the observed equivalent widths of the H/3, [OIIIjAA 4959, 5007 and Hct lines and 
continuum fluxes from free-free and free-bound processes, and convolving with the V and I 
bandpasses, they were able to reproduce the (V-I)o color range observed for the underlying 
component. We note here that these very blue colors arc not the results of some observational 
or reduction procedure artifact. We have obtained HST WFPC2 V and I images of the BCD 
Mrk 996 in exactly the same observational conflguration, and reduced them in exactly the same 
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manner (Thuan et al. 1996a). Yet, Mrk 996 shows a red underlying disk with a system of red 
globular clusters ((V-I)o ~ 0.9). 

3.3. A primordial neutral hydrogen intergalactic cloud? 

Here we present indirect evidence which also supports the youth hypothesis of SBS 0335-052. 
We use the BCD as a background light source shining through the HI envelope in which it is 
embedded (Thuan et al. 1996b) to probe the physical conditions of the surrounding neutral gas. 
The Lya line seen in absorption would give the column density of atomic hydrogen, while the 
01 A1302 absorption line would give the column density of the most abundant heavy element 
which remains neutral in the HI cloud. This would allow us to set a limit on the 0/H abundance 
ratio in the neutral gas. 

The spectrum showing the spectral region around the resonance line of oxygen at 1302 A 
is shown in Figure 6b. The 01 line is not seen in absorption. We can place an upper limit of ~ 
200 mA (corresponding to 2 times the rms noise) for its equivalent width, which corresponds 
to a column density of 2.7x10^"^ O atom cm~^, using the relation given by Kunth et al. (1994) 
between the equivalent width and the column density for an unsaturated line. Since the HI 
column density is ~ 10^^ H atom cm~^ from the VLA map (Thuan et al. 1996b), the oxygen to 
hydrogen abundance ratio is < 2.7x10"^. This is more than 3000 times smaller the N(0)/N(H) 
of Orion equal to 8x10"*^ and more than 3 times lower than the oxygen abundance measured 
in the HI envelope of IZwl8, where the 01 A1302 line was clearly detected in absorption ( 
Kunth et al. 1994). This observation provides further support for the youth hypothesis of 
SBS 0335-052. The HI cloud in which SBS 0335-052 is embedded appears to be composed of 
primordial pristine gas, uncontaminated by metal-enrichment from previous star formation. 

4. DISCUSSION 

Prom the (V-I)o colors of the super-star clusters and using Leitherer & Heckman (1995)' 
models, we have concluded that the present burst of star formation did not start until < 3-4 
Myr ago. From the very blue and filamentary structure of the underlying low-surface-brightncss 
component, we have concluded that the latter has not a stellar but gaseous origin. Because 
SBS 0335-052 does not appear to possess an underlying older stellar population, the present 
starburst is probably the first one in the galaxy's history, and the BCD is a good candidate for 
being a young galaxy. Is this conclusion consistent with other observational constraints? 

4.1. Dust 

Wc have shown that dust is clearly present in the BCD (see the color map in Figure 3). 
Can it be made in the short time scale implied by the young galaxy scenario? There is some 
observational evidence that dust formation can occur in a rather short time scale in supernova 
ejecta. Lucy et al. (1995) have attributed the brightening of SN 1987A in the IR after day ~ 
450 to condensation of small sificate grains. The dust condensation efficiency is rather low ( 
< 10~^ ), but an efficiency of ~1 is possible if the dust is clumped. Thus, while the physics 
of dust formation are still not very well known, the SN 1987A observations do not rule out 
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the possibility that enough sihcate grains can condense out of supernova ejecta to account for 
the dust observed in SBS 0335-052. The BCD has not been detected in any of the four IRAS 
bands. If we adopt as upper limit a flux of 0.5 Jy at 100 microns, then using the formula in 
Sauvage & Thuan (1994) and adopting their mean value f(100 microns)/ f(60 microns) ~ 2 for 
BCDs, we derive an upper hmit of ~ 1.6 x 10^ Mq for the warm dust mass. 

4.2. The supershell 

We next ask whether the size of the shell is compatible with a time scale of ~ 4 Myr. The 
radius of the shell is ~ 380 pc, which qualifies it to be a 'supershell' according to the definition 
of Tenorio-Tagle & Bodenheimer (1988). To put constraints on the age of the shell structure, 
we use the models of shell formation of McCray & Kafatos (1987). These take into account the 
effects of supernova as the encircled stellar clusters age. We need some estimate of the density 
of the ambient ISM. The HI mass of the whole cloud (which is much larger than that associated 
with the BCD) is ~ 10*^ M©, and the VLA map (Thuan et al. 1996b) shows that most of it is 
within a projected region 31 kpcx7 kpc in size (the whole cloud is larger, extending up to 67 
kpcx24 kpc, but the outer parts have a very low density and do not contribute significantly 
to the total mass). Assuming a cyhndrical shape for the HI cloud, an inchnation angle of 
degree and a uniform density gives a mean density of ~ 0.01 H atoms cm~^. As SBS 0335-052 
is associated with one of the two peaks in HI column density which is ~ 100 times higher than 
the mean, we obtain a density of ~ 1 H atom cm~^ for the ISM around the BCD. This is 
comparable to the density of ~ 1.5 H atoms cm~^ found by Hunter & Thronson (1995) for 
the ambient ISM of IZwl8. If we use the analytical expressions of McCray & Kafatos (1987) 
and assume that each supernova explosion produces an energy of 10^^ ergs, then to produce a 
supershell of radius ~ 380 pc in a time of ~ 4 Myr, then a total number of stars more massive 
than 7 Mq (all those massive enough to undergo a type II supernova) of ~ 14400 is required. 
This number is not unreasonable as the total H/3 luminosity (~ 4 xlO^° erg s^^) gives an equal 
number of ~ 10000 07 stars. The models predict a shell expansion velocity of ~ 56 km s~^. 

4.3. The star formation rate 

We can obtain an estimate of the SFR from the total Ho; luminosity. Using Melnick et al. 

(1992)' measurements, we obtain L(Ha) ~ 5.2x10^° erg s~^. This is a lower limit as Melnick et 
al. observed only SSCs 1, 2, 4 and 5, but not the fainter SSCs 3 and 6 which do not contribute 
significantly to the total H« fiux ( less than 10%). Extrapolating the Salpeter IMF down to 
stars of mass 0.1 Mq, the SFR is equal to 7.07x10"^^ L(Ha) M© yr'^ (Hunter & Gallagher 
1986). Since most of the ionizing photons are produced by 30-60 M© stars with hfetimes < 
3x10^ yr, this is the current instantaneous SFR. We obtain a SFR of ~ 0.4 Mq yr~^ for the 
BCD. This derived SFR is one order of magnitude larger than that in IZwl8 ( ~ 0.04 Mq yr^^ 
) , and is in the upper range of SFRs derived for BCDs (Fanelli, O'Connell & Thuan 1988). 
The total stellar mass formed during 4 Myr is thus ~ 1.6x10^ Mq. The total oxygen mass 
produced is ~ 6560 M© (Lequeux et al. 1981). 



4.4. The Lyman a emission 
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To obtain insight into the problem of the non-detection of primeval galaxies at high redshifts 
by Lya searches, we have obtained with HST a GHRS spectrum in the spectral region around 
Lycu at the redshift of SBS 0335-052. The spectrum is shown in Figure 6a. The most striking 
feature in this figure is the conspicuous absence of strong narrow Lya emission in SBS 0335- 
052, just as in IZwl8 (Kunth et al. 1994), and in contradiction with the predictions of models 
of young dust-free galaxies undergoing massive star formation (Chariot & Fall 1993). This 
problem is not new. Previous Lyo; studies of blue compact galaxies ( Deharveng et al. 1985, 
Hartmann et al. 1988, Terlevich et al. 1993, Giavalisco, Koratkar & Calzetti 1996) had already 
shown that Lya emission is very weak or absent in the UV spectra, despite the extremely 
strong optical emission lines. The favored explanation for the reduction of Lya fluxes from 
recombination values is absorption of multiply scattered Lya photons by dust in the HII region 
or in the surrounding HI envelope. This would imply increasing Lya fluxes in galaxies with 
decreasing metallicities, since presumably low-metallicity objects contain less dust, and hence 
suffer less destruction of Lya photons. The observational evidence is mixed. Terlevich et 
al. (1993) did find such a correlation, while Giavalisco et al. (1996) found none. The latter 
authors concluded that the ISM in BCDs is highly inhomogeneous and that the transport of 
Lya photons is primarily controlled by the ISM geometry rather than by the dust amount. The 
absence of strong narrow Lya emission in the two most metal-deficient BCDs known, IZwl8 
(Z0/5O) and SBS 0335-052 (Z0/4I), supports the lack of correlation between Lya fluxes and 
metallicities found by Giavalisco et al. (1996). 

In the case of SBS 0335-052, dust absorption may play some role in Lya photon destruction 
since dust is clearly present in the BCD (Figure 3). But it may not play the dominant role. Most 
of the light in the 2"x2" aperture comes from SSCl which, as we have argued, does not suffer 
from extinction. This conclusion is reenforced by the comparison of the observed continuum 
level near Lya (~ 7.5x10"^^ erg s~^ cm~^) with that near H/3 in the same 2"x2" aperture 
(Izotov et al. 1996). The ratio is ~ 100, close to the value predicted for young dust-free 
starbursts ( Chariot & Fall 1993). We favor therefore multiple scattering of Lya photons in an 
optically thick gas as the main mechanism of Lya attenuation. There is observational evidence 
that the gas surrounding SBS 0335-052 is optically thick for Lya photons. Figure 6a shows 
that, contrary to the situation in IZwl8 (Kunth et al. 1994) , the intensity of the central part of 
the Lya absorption caused by the HI envelope does not go to zero. In fact, superposed on the 
broader absorption, there is a broad (FWZI=20A) low-intensity emission with a total flux of 
2.6x 10~^^ erg cm~^ s~^. Broad emission is expected if resonant scattering occurs in a medium 
with optical depth of ~ 10^ for Lya photons (Bonilha et al. 1979). This optical depth would 
correspond to a HI column density of ~ 10^^ H atoms cm~^, which is about the observed value 
in the VLA map (Thuan et al. 1996b). Multiple scattering redistributes Lya emission over the 
whole extent of the HI cloud which is much larger than the aperture used. By using a very 
small aperture, most of the re-emitted photons are lost. This explains why the Lya-to-H/3 ratio 
in the BCD (~ 0.45) is so small compared to the case B recombination value of 33. 

The geometry of the HI cloud may also play a role. In SBS 0335-052, the HI cloud is seen 
nearly edge-on, as suggested by its flattened structure (Thuan et al. 1996b). Chariot & Fall 
(1993) have shown that, because photons escape more easily in a direction perpendicular to the 
plane than along it, the ratio of the Lya intensity to its mean intensity changes from a value 
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of 2.3 for face-on galaxies to 0.0 for edge-on galaxies. 

What are the implications of the above observations concerning the problem of the non- 
detection of high-redshift primeval galaxies in Lya searches? Dust is clearly present in a very 
metal-deficient environment such as in SBS 0335-052, so that dust absorption may play a 
role. Geometrical effects may also contribute. Galaxies seen face-on are easier to detect than 
those seen edge-on. Also we expect mctal-rich galaxies which have undergone repeated bursts 
of star formation to have a more eroded HI envelope with more holes through which Lya 
photons can escape. They would be more detectable than young very metal-deficient galaxies 
undergoing their very first burst of star formation with a less porous HI envelope. Lequeux et 
al. (1995) have discussed another mechanism which allows the escape of Lya photons. They 
detected attenuated Lya emission in the relatively metal-rich ( Z ~ Z0/3 ) BCD Haro 2. Some 
Lya photons could escape because they are moving at a different velocity from that of the 
surrounding HI gas, the velocity difference ( ~ 200 km s~^ ) being caused by a galactic wind. 
In summary, the escape and detection of Lya photons depends not only on the amount of dust 
present, but also on the geometry and velocity structure of the surrounding HI envelope. 

4.5. Comparison with other nearby young galaxy candidates 

Several other nearby objects have been proposed as possible young galaxies. This includes 
IZwl8, the only BCD known to be more metal-deficient than SBS 0335-052, and which also 
does not possess a low-surface-brightness underlying old stellar population (Kunth et al. 1988, 
Hunter & Thronson 1995). These two BCDs with no apparent underlying older stellar compo- 
nent constitute the exception rather the rule. 

The HI cloud associated with IZwlS is much smaller in intrinsic size ( 6 kpc x3 kpc, 
Viallefond, Lequeux & Comte 1987) than the one in which SBS 0335-052 is embedded. For 
both BCDs, the starburst region is associated with, but does not coincide exactly with the 
HI column density peak. This offset of ~ 1 kpc is a general feature of all BCDs with HI 
interferometric maps ( Viallefond & Thuan 1983). The HI cloud associated with SBS 0335-052 
has another HI peak at ~ 1.5 ' west (~ 24 kpc). There is a faint (V ~ 19 mag) compact 
(2.8"x2.2") dwarf galaxy at the location of the second HI peak. Its redshift ( v= 3990±40 km 

) is about the same as that of SBS 0335-052. Its spectrum shows Ha, H/3 and [OIII]A5007 
emission, but has lower excitation than that of the BCD ( Pustilnik et al. 1996). 

The very large size (64 kpc x 24 kpc) of the HI cloud of SBS 0335-052 is very unusual for 
BCDs. The typical size is more like a few kpc in each dimension (Viallefond & Thuan 1983). 
This large cloud is more reminiscent of the intergalactic cloud HI 1225-1-01 with a size of ~ 100 
kpc (adopting a distance of 10 Mpc, Salzer et al. 1991). Both clouds have HI masses of ~ 10^ 
Mq. There are also 2 HI peaks in HI 1225-1-01 separated by ~46 kpc, one of which is associated 
with a faint ( ~ -14.0 ) and metal-deficient ( Z ~ Z0/I8 ) dwarf irregular galaxy. As for 
SBS 0335-052, the galaxy is very blue ( B-V= 0.10, U-B= -0.57 ) and most ( > 72 percent ) 
of the B light comes from stars with ages < 40 Myr. Thus, although the burst in SBS 0335-052 
is younger than in HI 1225-1-01 , their HI envelopes have remarkably similar properties. It is 
not yet clear why star formation has been stunted in these HI clouds for so long. In the case 
of SBS 0335-052, a mean density of ~ 0.01 H atom cm~^ leads to a collapse time of only ~ 
5x10^ yr, much shorter than the Hubble time. 
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5. SUMMARY 

We have analyzed WFPC2 V and I images and GHRS UV spectra obtained with the refur- 
bished HST of the extremely metal-deficient ( Z = Z0/4I ) blue compact dwarf (BCD) galaxy 
SBS 0335-052. We have obtained the following results: 

1) The underlying extended low-surface-brightness component of the BCD is extremely 
blue (-0.34 < (V-l)o < 0.16 ). Its color profile is nearly flat. If stars are responsible for this 
low-surface-brightness emission, they cannot be older than ~ 10^ yr. However, the irregular 
and filamentary structure of the underlying component supports the case for gaseous emission 
and no underlying older stellar population. The very blue colors are consistent with gaseous 
emission. Thus SBS 0335-052 is very likely a nearby young galaxy in the sense that it is 
undergoing its very first burst of star formation. This conclusion supports the prediction of 
Izotov et al. (1990) about the youth of SBS 0335-052 on the basis of its very low metaUicity. 

2) The 01 A1302 line is not detected in absorption in the GHRS spectrum, setting an upper 
limit for N(0)/N(H) in the HI envelope of ~ 3xl0~^, more than 3000 times smaller than the 
value in Orion and more than 3 times lower than the oxygen abundance measured in the HI 
envelope of IZwl8, and consistent with the youth hypothesis of SBS 0335-052. 

3) All star formation occurs in six very blue super-star clusters (SSC) with ages < 3-4 
Myr. SSC formation appears to be related to starburst events in some, but not all, BCDs. The 
triggering of star formation in SBS 0335-052 is not due to tidal interactions, but is probably 
caused by stochastic collisions of gas clouds within the large HI envelope. 

4) Dust is clearly present and is spatially mixed with the SSCs. While the brightest SSC 
does not suffer from extinction, there is a systematic increase in reddening away from it. Such 
dust can probably condense out in supernova remnants in a time scale < 4 Myr. 

5) There is a supershell of radius ~ 380 pc delineating a large supernova cavity. Its observed 
properties can be accounted for by models of shell formation with a time scale of ~ 4 Myr. 

6) The star formation rate (SFR) is ~ 0.4 Mq yr~^, in the upper range of SFRs observed 
in BCDs. 

7) Strong narrow Lycu emission is not detected. There is instead low-intensity broad (FWZI= 
20A) Lya emission superposed on broader Lya absortion from the HI envelope. The low- 
intensity broad emission is caused by resonant scattering in a dense ambient medium, with 
optical depth of ~ 10^ for Lya photons, corresponding to the observed HI column density of ~ 
10^^ H atom cm-^ 

While dust may play some role in the destruction of Lya photons, multiple scattering which 
removes Lya photons from the small HST aperture and redistributes them over the whole HI 
cloud, and geometrical effects ( the HI cloud is seen edge-on) appear to be the main mechanisms 
for removal of Lya photons in SBS 0335-052. Lya detection of high-redshift primeval galaxies 
depends thus not only on the amount of dust, but also on the geometry and velocity structure 
of the surrounding HI envelope. 
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FIGURE CAPTIONS 

Fig.l. (a) A 1800 s V WFPC2 image of the blue compact dwarf galaxy SBS 0335-052 
taken with the refurbished HST, with the contrast adjusted to show the high surface brightness 
super-star clusters (SSCs). At a distance of 54.3 Mpc, 1" corresponds to a linear size of 263 pc. 
The orientation is the same as in Figure lb. (b) A V contour map showing the labeling of the 
six SSCs listed in Table 2. A supershell delineating a large supernova cavity of radius ~ 380 
pc is clearly seen at the NE end of the galaxy. 

Fig. 2. The (V-I)o color corrected for foreground extinction of the six super-star clusters 
(SSCs) listed in Table 2 as a function of their separation from SSCl. There is a systematic 
increase in reddening with increasing distance from SSCl, along the SE-NW direction. 

Fig. 3. (V-I) grey scale color map. Blue is dark and red is light. The supershell is very blue 
while the red dust patches are seen to increase in importance away from super-star cluster 1 in 
the SEl-NW direction. Two background red galaxies can be seen to the south of SBS 0335-052. 

Fig. 4. The same V image as in Figure la but with the contrast adjusted to show the 
low-surface-brightness underlying component of SBS 0335-052. The latter has a blotchy and 
irregular appearance, suggesting that the light is not emitted by stars but by ionized gas. The 
supershell delineating the large supernova cavity is clearly seen. 

Fig. 5. (a) and (b) V and / surface brightness profiles as a function of log r where r is the 
equivalent radius. A power law is a straight line in this coordinate system. The error bars 
take into account photon statistics. The large error bars in the radius range -0.45 < log r 
< -0.05 are caused by light contamination from the shell structure, (c) (V-I) color profile. 
It is approximately fiat and shows an extraordinary blue color for the underlying low-surface- 
brightness component. 

Fig. 6. GHRS spectra obtained through a 2"x2" rectangular aperture with HST and centered 
on super-star cluster 1. The aperture includes also super-star clusters 2 and 3. The spectral 
resolution is 0.12 A. (a) shows the spectral region around Lya. The strong geocoronal Lya line 
can be seen at 1216 A. The redshifted location of the Lya fine in SBS 0335-052 (v = 4043 km 
s~^) is marked. There is low-intensity broad (from 1220A to 1240A) hya emission superposed 
on even broader absorption. The continuum level is 7.5x10"^^ erg cm~^ s~^ A~^. (b) shows 
the spectral region around the 01 A1302 line. Its redshifted location is marked. The line is not 
detected and an upper hmit of - 200 mA can be set for its equivalent width. 
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Table 1. Observational characteristics of SBS 0335-052 



Parameter Value 



a(1950) 03''35""15.2^ 

5(1950) -05°12'26" 

F, b^^ 193°, -45° 

d25" 12.6" 

D (Mpc)'' 54.3 

d (kpc) 3.29 

V" 16.65 ±0.01 

V-I" -0.23±0.01 

Mb -16.70 

whi (km s-^Y 4043 ±10 

Av2o (km s"^)^ 105 ±14 

Avso (km s-^)= 83 ± 9 

¥hi (10^ M0 Mpc-^)'^ 0.34 

Mhi (10^ MqY 0.99 

Lb (109 Lq)'^ 1.1 

Mbj/Lb (Mq/Lq) OJO 



" from this paper. d25 is derived from Figure 2. V and I measured in a 12" diameter aperture 
centered on the brightest knot. 

redshift distance corresponding to Hq = 75 km s~^Mpc~^ and a systematic velocity of 4076 
km s~^, corrected to the Local Group velocity centroid following Thuan & Martin (1981). 

HI measurements obtained with the Nangay radio telescope by Thuan et al. (1996); \hi 
is the heliocentric velocity. 

Adopting MB(sun) = 5.48 mag. 
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Table 2. Magnitudes and colors of stellar clusters 



No VI /g {V - 1)1 



1 18.93± 0.01 19.54± 0.01 -0.61± 0.02 

2 19.24± 0.02 19.63± 0.02 -0.39± 0.02 

3 20.55± 0.08 20.81± 0.10 -0.26± 0.12 

4 19.71± 0.02 19.86± 0.02 ~0.15± 0.03 

5 19.50± 0.03 19.36± 0.03 0.14± 0.04 

6 21.96± 0.11 21.65± 0.09 0.31± 0.14 



"corrected for foreground extinction (Ay^O.lO mag, A/= 
0.06 mag, Burstein & Heiles 1982; Whitford 1958). 
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